Possibility of superconductivity in graphite intercalated with alkaline earths 
investigated with density functional theory. 

Matteo Calandra and Francesco Mauri 
Institut de Mineralogie et de Physique des Milieux condenses, 
case 115, 4 plcice Jussieu, 75252, Paris cedex 05, France 
(Dated: February 6, 2008) 

Using density functional theory we investigate the occurrence of superconductivity in ACe with 
A=Mg,Ca,Sr,Ba. We predict that at zero pressure, Ba and Sr should be superconducting with 
critical temperatures (Tc) 0.2 K and 3.0 K, respectively. We study the pressure dependence of Tc 
assuming the same symmetry for the crystal structures at zero and finite pressures. We find that 
the SrCe and BaCe critical temperatures should be substantially enhanced by pressure. On the 
contrary, for CaCe we find that in the to 5 GPa region, Tc weakly increases with pressure. The 
increase is much smaller than what shown in several recent experiments. Thus we suggest that in 
CaCe a continous phase transformation, such as an increase in staging, occurs at finite pressure. 
Finally we argue that, although MgCe is unstable, the synthesis of intercalated systems of the kind 
MgxCai-xCy could lead to higher critical temperatures. 

PACS numbers: 63.20.Kr, 63.20.Dj , 78.30.Er, 74.70.Ad 



I. INTRODUCTION 

It has long been known that Graphite intercalated 
compounds (GICs) can display a superconducting be- 
havior at low temperature^. However, until the dis- 
covery of ytterbium and calcium intercalated graphiteSi^ 
(Tc(YbC6)=6.5 K and Tc(CaC6)=11.5 K), the critical 
temperatures achieved were typically inferior to 5 Kelvin. 
Very recently, it has been shown that even higher critical 
temperatures (up to 15.1 K) can be achieved in CaCe 
applying hydrostatic pressure (up to 8 GPa)^. This is 
presently the highest Tc reported in a GIC. 

From a practical point of view GICs are appealing 
since carbon is a very versatile element. Beside the fact 
that the three most known carbon crystalline structures 
(Fullerene, Graphite and Diamond) display supercon- 
ducting behavior upon intercalation^'^ or doping^, car- 
bon is also interesting for the possibility of making nan- 
otubes which can display a superconducting behavior^-. 
Moreover the number of reagents which can be inter- 
calated in graphite, by using chemical methods or high 
pressure synthesis, is very large. Consequently finding 
GICs with even higher critical temperatures is not a re- 
mote possibility. Therefore theoretical calculations can 
be a precious tool for the investigation of superconduc- 
tivity in GICs. 

The nature of superconductivity in intercalated GICs 
is still controversial. Due to their layered nature and 
to the presence of an intercalant band at the Fermi level 
(Ej), Csanyi et al.^ suggested a non-conventional exciton 
pairing mechanism^. On the contrary MaziniS proposed 
an electron-phonon mechanism mainly sustained by Ca 
vibration. In a previous work^^ we suggested that the 
pairing is mediated by the electron-phonon interaction. 
In particular the carriers are mostly electrons in the Ca 
Fermi surfaces coupled with Ca in-plane and C-out of 
plane vibrations. The calculation of the isotope effect 
coefficient showed that the contribution of Ca in-plane 



vibrations and C out-of-plane vibrations to superconduc- 
tivity is comparable. 

Experimental data seem to confirm that the pairing 
in CaCe and YbCe is indeed due to the electron-phonon 
interaction, but several open questions remain. From 
the exponential behavior of the penetration depth in 
CaCe Lamura et al}'^ deduced an isotropic gap of mag- 
nitude A(0)=1.79± 0.08 meV. A similar isotropic gap 
(A(0)=1.6± 0.2 meV) has been measured by scanning 
tunneling spectroscopy^^. The corresponding values of 
2A(0)/Tc are in agreement with the BCS theory. Simi- 
lar conclusions have been obtained from the specific heat 
jump at superconducting transition^^. Thermal conduc- 
tivity data in the presence of a magnetic field^^ indicate 
that in YbCe the order parameter has s-wave symmetry 
too and exclude the occurrence of multiple gaps. Iso- 
tope effect measurements^^ show a huge Ca isotope co- 
efficient, a(Ca) = 0.5, in disagreement with theoretical 
calculations^^. This is surprising, since the gap and spe- 
cific heat data are correctly described by DFT calcula- 
tions, meaning that the calculated electron-phonon cou- 
pling is probably correct. Even more interesting is that 
the total isotope effect would be probably larger than 0.5, 
although C-isotope effect measurements are necessary to 
confirm this and to judge the validity of the measure- 
ments of a(Ca). In any case, the large Ca-isotope coef- 
ficient, the measured superconducting gap and the jump 
of the specific heat at the transition they all go in favor of 
a phonon mediated mechanism with, most likely, a single 
s-wave gap. 

In this work we want to push a step forward the pre- 
diction that can be made by the electron-phonon theory. 
This is important since it may allow to identify GICs 
with higher critical temperatures and it also represents 
a significative benchmark for DFT simulations. It was 
noted^ that all the superconducting GICs possess an in- 
tercalant Fermi surfaces at Ej. This fact is relevant for 
both pairing mechanisms proposed. In the case of a con- 
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ventional electron-phonon mechanism the electrons in the 
intercalant Fermi surface are the ones strongly coupled 
to the phononsil. Thus we study the possible occur- 
rence of superconductivity in graphite intercalated with 
alkaline earths (ACe with A=Ba,Sr,Mg). Ah these GIGS 
have an intercalant Fermi surface at Ej so they are good 
candidates for superconductivity. We predict the critical 
temperatures for BaGe and SrGe in the framework of the 
electron-phonon coupling theory. 

As mentioned before, the critical temperature of GaGe 
and YbGe is enhanced with pressure. Resistivity mea- 
surements under pressure show that, at 8 GPa, GaGe 
undergoes a structural phase transition to a new super- 
conducting phase with a lower critical temperature. The 
new structure seems to be stable at least up to 16 GPa. 
In other successive works magnetic^^'^^ and resistive^^ 
measurements were performed in a much smaller range 
of pressures (0-1.6GPa) and the behavior observed was 
confirmed. Pressure also increases the critical tempera- 
ture of YbGe up to approximatively 7.0 K at ^2.0 GPa. 
In this case too a structural transition is seen towards a 
new superconducting phase with lower Tc. In both YbGe 
and GaGe the dependence of Tc is linear with similar 
ATc/AP (0.4 for YbGe and 0.5 for GaGe). The fact that 
Tc can be enhanced with pressure suggests that this can 
be a general mechanism for superconducting in GIGs. It 
is then important to study superconductivity as a func- 
tion of pressure in BaG6,SrG6 and finally GaGe- 

After illustrating in section ^ and IIIII the technical 
parameters and the lattice structures used in the sim- 
ulations, we study the electronic structure (sec. IIVI ), 
the phonon dispersions (sec. and the superconduct- 
ing properties (sec. IVI|) of Alkaline-earth intercalated 
graphite at ambient pressure and at finite pressure. Par- 
ticular emphasis is given to the case of GaGe under pres- 
sure (sec. IVII|) . 



II. TECHNICAL DETAILS 

Density Functional Theory (DFT) calculations are 
performed using the espresso code^^ and the general- 
ized gradient approximation (GGA)^^. We use ultrasoft 
pseudopotentials^ with valence configurations 3s^3p^4s^ 
for Ga, 4s24p^4d^5s^5p0 for Sr, Ss^Sp^Sd^s^Gp^ for Ba, 
and 2s^2p^ for G. For Mg we use Troullier-Martin^^ pseu- 
dopotentials with configuration 3s^-^3p^3d^. The wave- 
functions and the charge density are expanded using a 
35 Ry and a 600 Ry cutoff. The dynamical matrices 
and the electron-phonon coupling are calculated using 
Density Functional Perturbation Theory in the linear 
response For the electronic integration in the phonon 
calculation (structure R3m) we use a =6x6x6 and 
Nk = 8x8x8 uniform k-point meshes and Hermite- 
Gaussian smearing ranging from 0.1 to 0.05 Ry. In or- 
der to obtain very accurate phonon-frequencies for the 
low energy modes (below 15 meV) it is crucial to use 
a large cutoff for the charge density (600 Ry at least) 



and a very high convergence threshold in the phonon 
calculations. For the evaluation of the electron-phonon 
coupling and of the electronic density of states we use 
Nk = 25 X 25 X 25 and A^"^ = 20 x 20 x 20 meshes respec- 
tively. For the A average over the phonon momentum q 
we use aA/'g = 4x4x4 q— points mesh. The phonon 
dispersion is obtained by Fourier interpolation of the dy- 
namical matrices computed on the Nq mesh. 



III. CRYSTAL STRUCTURE 

All the considered compounds are stage 1^ at zero pres- 
sure. The atomic structure^ of GaGe involves a stacked 
arrangement of graphene sheets (stacking AAA) with Ga 
atoms occupying interlayer sites above the centers of the 
hexagons (stacking aPj). The crystallographic struc- 
ture is R3m where the Ga atoms occupy the la Wyckoff 
position (0,0,0) and the G atoms the 6g positions (x,- 
x,l/2) with x= 1/6. SrGe and BaGe have a slightly 
different crystal structure^^ involving a stacked arrange- 
ment of graphene sheets (stacking AAA) with Sr and 
Ba atoms occupying interlayer sites above the centers of 
the hexagons with an a/3 stacking. The crystallographic 
structure is P63/mmc where the Sr,Ba atoms occupy the 
2c Wyckoff position (1/3,2/3,1/4) and (1/3,2/3,3/4) and 
the G atoms the 12i positions (1/3,0,0). The experimen- 
tal in plane lattice parameter a and the interlayer spacing 
c for the three structures are illustrated in tabled 

Material Stacking a,c a,c (LDA) a,c (GGA) 

BaCe AaAp 4.302 5.25 4.280 5.00 4.350 5.20 

SrCe AaAp 4.315 4.95 4.285 4.80 4.325 5.00 

CaCe AaA;^A7 4.333 4.524 4.29 4.36 4.333 4.51 

MgCe 4.317 3.75 4.35 3.95 



TABLE I: Experimental (first two columns) structural pa- 
rameters (Angstrom) of BaCe, SrCe, CaCe- With c we indi- 
cate the interlayer spacing, namely the distance between two 
graphene layers, while a is the in plane lattice parameter. In 
the last two columns we report the parameters of the theo- 
retical (LDA and GGA) R3m structure (stacking AaApA^) 
having the same a and the same interlayer distance of the 
experimental structure. This rhombohedral structure is con- 
sidered in all the calculations. Since MgCe has has never been 
synthesized we only report its theoretical minimized structure 
with symmetry R3m 

Even if the structural symmetry of SrGe and BaGe are 
different from the one of GaGe, in this work we consider 
the same rhombohedral symmetry group (namely R3m 
with stacking AaA^Aj ) for all the considered GIGs. We 
do not expect this assumption to affect qualitatively the 
calculated electronic and phonon properties since the two 
structures differ only for large distance neighbors. Indeed 
the distances of the first and second nearest neighbors are 
the same. The differences between the metal lattice sites 
in the AaA/SAj and in the AaA/3 structures is equivalent 
to those existing in the fee and in the hep structures. 
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System 


Pressure (GPa) 


a 


c 


BaCe 


8 


4.287 


4.925 


SrCe 


8 


4.295 


4.750 




16.5 


4.265 


4.524 


CaCe 


3 


4.333 


4.524 




5 


4.317 


4.349 




6 


4.314 


4.330 




7 


4.310 


4.301 




8 


4.307 


4.290 




9 


4.303 


4.270 




10 


4.300 


4.250 



TABLE II: Theoretical (GGA) Structure of Alkali- earth GICs 
under pressure 

In the last three columns of tabled we report the theo- 
retically minimized parameters assuming the R3m struc- 
ture and using the local density approximation (LDA) 
or the Generalized- Gradient approximation (GGA). The 
interlayer distance between graphite layers (c) and the 
in-plane lattice parameter (a) calculated with GGA are 
in very good agreement with experiments. The equilib- 
rium LDA a and c parameters are slightly compressed 
respect to the CGA and experimental values. 

Since MgCe has never been synthesized, we also as- 
sume that it crystallizes in the R3m structure. The the- 
oretical equilibrium parameters are reported in table ffl 
We have verified that this structure is unstable since its 
energy in CGA is lower by 0.016 Ry/(Cell CaCe) then 
the ones of Magnesium and Graphite separated. 

In Alkali-earths GICs as the atomic number of the in- 
ter calant (Z) is reduced the c parameter is also reduced. 
To disentangle the two effects we consider BaCe, SrCe 
and CaCe under isotropic pressure. Since the structure 
of these systems at a given pressure is not known, we 
minimize the R3m structure at a given isotropic pres- 
sure. The lattice parameters obtained are illustrated in 
tabled It should be noticed that for SrCe the pressure 
of 16.5 GPa has been chosen since the c parameter has 
the same value as in CaCe at zero pressure. 

When possible, we always used the experimental a 
and c parameters in the electronic structure calculations. 
When lacking, we used the GGA minimized parameters. 



IV. ELECTRONIC STRUCTURE 

The zero-pressure DFT-band-structures of the 
Ba,Sr,Ca,Mg intercalated compounds are illustrated 
in fig. n All the considered alkali-earths intercalated 
GICs have at least one intercalant band at the Fermi 
level. This was proposed to be^'^^ a necessary condition 
in order to have superconductivity with a reasonable 
critical temperature in GICs. 

The considered GICs are layered structures, but their 
band structure are clearly three-dimensional. This was 
understood in the case of CaCe^^ since the rhombohedral 



angle (49.55*^) is not too far from the one corresponding 
to the fee structure (60^). The rhombohedral angles for 
BaCe, SrCe and MgCe are 43.47^45.81^ 55.37^ respec- 
tively. Thus the rhombohedral angle is larger as c is 
reduced. This is even more evident comparing the be- 
havior of the band-structure as a function of pressure in 
fig. 121 for SrCe- Note in particular, that at 16.5 GPa 
(where the c parameter is the same as that of CaCe at 
zero pressure) the band structure is very similar to the 
one of CaCe at zero pressure. The larger bandwidths of 
the intercalant band as c is reduced reflects the higher 
three dimensional character of the structure. 

The electronic density of states for the four considered 
compounds are very similar (see fig. ^ and also their 
values at the Fermi level are very close. No qualitative 
difference in the superconducting behavior of these sys- 
tems can be attributed to the number of carriers at the 
Fermi energy (i.e. to the different value of N(0)). 

The atomic-projected density of states is cal- 
culated using the Lowdin populations, Pr]{e) = 
7^ Ekn K</>^IV^kn)p(^(ekn " e)- In this expression |(/)^) = 
5^77/ [^~^^^]^,^' 1^^') orthonormalized Lowdin or- 

bitals, are the atomic wavefunctions and Sn^rj' = 

. The density of states projected over the inter- 
calant states has very similar value at Ej for the three 
case and it is slightly reduced with decreasing Z. No 
particular differences respect to the case of CaCe are 
observed^^. In fig. Q (right) the size of the dots rep- 
resents the projection over the intercalant atomic states. 

The three Fermi surface sheets of BaCe, SrCe and 
CaCe are shown in fig. |2l The sheets have been identified 
by the corresponding band index. As c is reduced pass- 
ing from BaCe to CaCe the main difference is the change 
in the third sheet of the Fermi surface which is indeed 
a warped cylinder in BaCe and it becomes spherical in 
CaCe- This is best seen analyzing the band structure in 
the TL direction (namely the kz direction). In BaCe and 
SrCe there are no bands crossing Ej along FL, while in 
CaCe the Ca-originated band crosses Ej. 

V. PHONON DISPERSION 

The phonon dispersion of all the considered GICs can 
be divided in three regions. The high energy region is 
composed by in-plane Carbon- vibration ( Cxy ), the in- 
termediate energy region by out-of-plane Carbon vibra- 
tion ( C^ ), and the low energy region is mostly composed 
intercalant vibration ( l^y for in-plane and I^ for out-of- 
plane, with I=Ba,Sr,Ca,Mg), as it was shown for CaCe^. 
Since the high energy branches {uc^u > lOOmeV ) are 
weakly affected by the type of intercalant, for the sake 
of clarity, in some case we only plot branches in the low 
energy region of the spectra. These branches are those 
having the greatest contribution to the electron-phonon 
coupling. 

The zero-pressure phonon dispersion of BaCe, SrCe 
and CaCe are illustrated in figs. EJandO As Z is reduced. 




FIG. 1: (Color online) Band structure (left) and density of states projected over atomic orbitals (right) for several Alkali- 
earths GICs. The zero energy corresponds to the Fermi energy. The size of the dots represents the percentage of intercalant 
component. As a reference, in CaCe the dot at ~ 0.6eV at the L point represents 95%. 



and consequently c is reduced, the lowest energy branch 
in the FX direction (which corresponds to in-plane mo- 
menta with kz = 0) is softened. This is evident in CaCe 
where an incipient Kohn anomaly is found at X. While 
the reduction of c or Z softens the Ixy, the intercalant 
out-of-plane vibrations are hardened. 

This effect is also apparent from the analysis of the 
phonon dispersion of SrCe and CaCe under pressure. In 
SrCe the lowest branch at 16 GPa in the FX direction 
is almost imaginary, meaning that at this pressure the 
system is close to a structural instability^^. This insta- 
bility is driven by the reduction in the c-lattice spacing. 
In MgCe, at the theoretical lattice parameters U we find 
imaginary phonon frequencies at X, so that the system 
is probably unstable. As a consequence we did not cal- 
culate the full phonon dispersion. 



pling Aqj, for a phonon mode u with momentum q: 



k,n,m 



l5'kn,k+qm I^'^(ekn)i5(ek+qm) 



(1) 

where the sum is over the Brillouin Zone. The matrix 
element is g^^ 

,k+qm ~ (kn| JT^/(^iiqzy |k + qm)/y^2aJqz7, 
where Uc^u is the amplitude of the displacement of the 
phonon and V is the Kohn-Sham potential. The aver- 
age electron-phonon coupling is A = X]q^ ^c^v/Nq and its 
value for the different compounds and for different pres- 
sures is given in table II VI 
The Eliashberg function 



1 



(^^F{U) = ^ Aq^CJq^^(cJ 



(2) 



VI. SUPERCONDUCTING PROPERTIES 

The superconducting properties of Alkali-Earths GICs 
can be understood calculating the electron-phonon cou- 



and the integral A(cj) = 2 duo' a^F{(ji}') /uo' are shown 
in fig. 13 for BaCe and SrCe- As for CaCe the largest 
contribution to A comes from modes below 75 meV. 

An estimate of the different contributions of the in- 
plane and out-of-plane vibrations of the different atoms 
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L r 



P=16.5 Gpa 

FIG. 2: Band structure of SrCe at different pressures. The 
zero of the energy corresponds to the Fermi energy 




FIG. 3: (Color onhne) Fermi surface sheets of BaCe, SrCe 
and CaCe at zero pressure. The sheets are plotted in order of 
increasing band index. 



to A can be obtained from the relation 



(3) 



q iajf3 




(r|,r|,r|) 



(T1-T1,0) 



(T1,0,0) 



where i^a indexes indicate the displacement in the 
Cartesian direction a of the i^^ atom, [Gci]iaji3 = 
Ek,n,m4^*a^j73^(ekn)^^(ek+qm)/[A^(0)A^/c], and = 
(kn|(5y/feq^Q, |k + qm)/A/2. The Cq matrix is the 

Fourier transform of the force constant matrix (the weakly affected by Z or c-axis reduction. 



FIG. 4: Phonon dispersions of BaCe, SrCe and CaCe at 
GPa. 



derivative of the forces respect to the atomic displace- 
ments). 

Using eq. 01 we decompose A restricting the summa- 
tion over z, a and that over on two sets of atoms and 
Cartesian directions. The sets are Cxy^ C^, Ixy^ and I^, 
where I={Ba,Sr,Ca}. The resulting A matrix for the dif- 
ferent GICs and as a function of pressure are given in 

table mn 

Except for of SrCe at 16.5 GPa, the off-diagonal matrix 
elements are negligible. In all the case the largest con- 
tributions to A comes from Cz and Ixy phonon modes. 
The coupling to these modes are enhanced as Z and c 
are reduced. A similar enhancement of Cz and Ixy oc- 
curs under pressure. When the pressure is too high, as 
in SrCe at P=16.5 GPa, the off-diagonal matrix elements 
increase and the attribution to Cz or to Ixy becomes ill- 
defined. Note also that the Cxy and the Iz modes are 



(ri,ri,ri) 



(T1-T1,0) 



(ri,ri,0) 



(T1,0,0) 



FIG. 6: Phonon dispersions of SrCe at 8 and 16.5 GPa. 



FIG. 5: (Color online) Low energy region of the phonon dis- 
persions of BaCe, SrCe and CaCe at GPa. The size of 
the blue dots (red squares) represents the amount of Ixy (I^), 
where I = Ba, Sr, Ca. vibrations. 



The critical superconducting temperature is estimated 
using the McMillan formula^^: 



M 

1.2 



1.04(1 + A) 
" X-^x* (1 + 0. 62A) 



(4) 



where /i* is the screened Coulomb pseudopotential and 

(uj) = ex/o^°°"'^MiogM/^^^ (5) 

the phonon frequencies logarithmic average. Results for 
(uj) and for Tc are shown in table|^using /i* = 0.14. For 
BaCe and SrCe, Tc increases as the interlayer spacing is 
decreased. For SrCe at 16.5 GPa the the result of McMil- 
lan formula is not correct. Indeed close to the transition 
the increase in A is only given by the softening of the 
Sixy vibration close to X. In McMillan formula the limit 
^qi/ implies Tc ^ as {u) since Aq^^ ^ oo as l/uj'^^. 
Unfortunately the use of McMillan formula in this limit is 
not correct since it deviates from the Migdal-Eliashberg 
results(ME), as known from the papers of McMillanS^ 
and of Allen and Dynes^^(see fig. 1 in^^). For this rea- 
son the calculation of Tc for SrCe at P=16.5 GPa should 
be considered not reliable (and as a consequence it is not 
given in table IIV|) . 



VII. CaCe UNDER PRESSURE 

The study of CaCe under pressure needs particular 
emphasis due to the interest motivated by the work of 
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FIG. 7: (a) Zero-pressure Eliashberg function, a^F{co), (con- 
tinuous line) and integrated coupling, X{co) (dashed) for BaCe, 
SrCe and CaCe- 



Gauzzi et alA showing a considerable increase of Tc 
which reaches its maximum of 15.1 K at 7-8 GPa. The 
increase of Tc has been confirmed by other two succes- 
sive works measuring Tc in a much smaller range of 
pressure^^'^^. At 8GPa the system seems to undergo 
a phase transition towards a different structure with a 
lower superconducting Tc- 
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To study the behavior of Tc under pressure we calcu- 
late the phonon dispersion and the electron-phonon cou- 
pling at 5GPa, for which a Tc = 14:K was detected in 
experiments. The comparison between the phonon dis- 
persions at OGPa and 5GPa is shown in fig. |H1 As can be 
seen, while the lowest Ca phonon mode is slightly soft- 
ened (mostly at X), the others Ca and Cz phonon modes 
are hardened. However the overall change is fairly small 
and the Eliashberg functions a'^F{uj) at and 5GPa have 
no sizable differences (see fig. EJ. At higher pressures the 
phonon frequencies become imaginary, similarly to what 
happens in SrCe at 16.5GPa. In particular the largest 
softening occurs at X for the lowest Ca mode and it 
becomes negative in the range of 7-10 GPa. The low- 
est phonon frequency at X is extremely dependent on 
the Hermite-Gaussian smearing used in the calculation 
and as a consequence much larger k-point mesh must 
be used to correctly identify the phase transition (using 
a 10X10X10 mesh and smearing 0.04 still gives results 
which are not converged). However it is clear that at a 
sufficient high pressure the system will become unstable 
and the lowest Ca phonon frequency at X will approach 
zero. 

More insight on the superconducting properties can 
be gained by computing a;iog and A at and 5 GPa. In 
this range of pressures these parameters have very similar 
values and consequently the critical temperature changes 
only slightly (see table HVj), from 11.03 K at OGPa to 
11.40 K at 5GPa. 

The calculation for Tc under pressure are in stark 
disagreement with a recent theoretical calculationii ob- 
tained using the same method (DFT) and the same code. 

The disagreement with the calculation in ref.^^ can be 
explained by the following. In refi the electron-phonon 
coupling is calculated close to the structural transition 
(at lOGPa) where the Ca phonon frequency strongly soft- 



ens. Just before the structural transition, c^ca ^ and 
consequently Aca ^ oo. In this limit (close to the tran- 
sition), t/ie behavior is highly non linear. On the contrary 
in ref.^^ the behavior is assumed linear between OGPa 
and lOGPa^^. Our calculation at 5.0 GPa giving essen- 
tially the same Tc as at GPa shows that this is not the 
case. 

We remark that in the calculation of the electron- 
phonon coupling at and 5GPa it is crucial to use at 
least a 4x4x4 q-point mesh. If the smaller 2x2x2 mesh 
is used, the electron-phonon coupling of the Ca modes is 
overestimated, since the point at X (which is included 
in the mesh and has a considerable electron-phonon cou- 
pling) has a too large weight. The consequent increase of 
Tc under pressure is much larger than in the 4 x 4 x 4x 
mesh. 

From the preceding discussion it follows that DFT cal- 
culations give a non linear behavior of Tc versus pressure, 
in stark disagreement with experiments. In particular 
the increase of Tc as a function of pressure is too weak 
when compared with experiments. In the next section 
we discuss what can be the origin of such disagreement. 



VIII. CONCLUSIONS 

In this work we theoretically investigated the occur- 
rence of superconductivity in graphite intercalated with 
Alkaline earths (in particular Ba,Sr,Ca,Mg). Most of 
these systems have been already synthesized (BaC6,SrC6 
and CaCe), while MgCe has been proposed^^ as a good 
candidate for superconductivity due to its light mass and 
small force constants, possibly leading to large electron- 
phonon coupling. 

We predict the critical temperatures of BaCe and SrCe 
to be 0.23 K and 3.02 K respectively. Moreover we also 
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FIG. 9: (Color online) Comparison between a^F(cj) and A(a;) at GPa and 5 GPa for CaCe 



predict a substantial increase in the critical temperature 
under pressure for these two systems, namely at 8 GPa 
the critical estimated temperatures for BaCe and SrCe 
are 0.81 K and 5.16 K respectively. We hope that these 
predictions can stimulate additional experimental work 
so that the nature of the pairing in GICs can be further 
elucidated. Moreover it would be important to judge 
the reliability of DFT in calculating the superconduct- 
ing properties of different GICs. Indeed, while supercon- 
ducting gap and specific heat measurements are in very 
good agreement with DFT calculations the only available 
measure of the Ca isotope effectii^ in CaCe is in stark dis- 
agreement with DFT predictionsii. The claim that the 
isotope effect coefficient of Ca is 0.5, if verified by the cor- 
responding measurement of the C isotope effect, would 
open new perspectives. 

We have also shown that MgCe is energetically unsta- 
ble against phase separation in Mg and graphite. More- 
over, even assuming that it can be synthesized using pres- 
sure methods, the minimized structure at zero pressure 
obtained using DFT has imaginary phonon frequencies 
at the point X. Nonetheless it is interesting to note that 
the softening at the point X occurs when the interlayer 
spacing is reduced in CaCe and is connected to a large 
electron-phonon coupling in this point. This information 
suggests that the Tc of CaCe can probably be raised syn- 
thesizing Mga^Cai-a^Cjy alloys. 

Concerning CaCe, a puzzling problem is the depen- 
dence of the critical superconducting temperature upon 
pressure. Experimentally Tc increases as a function of 
pressure, but as we have shown in this work, Tc in- 
creases much faster than what DFT calculations predict. 
Further work is necessary to explain the origin of this 



disagreement, mainly on the experimental side. Indeed 
diffraction data as a function of pressure are absolutely 
necessary. If there is a phase transformation upon pres- 
sure, the structure minimized with DFT at finite pres- 
sure would not be correct. The fact that the Tc ver- 
sus pressure curves are very smooth seems to exclude an 
abrupt transition. A possibility is that a staging transi- 
tion occurs in CaCe under pressure. This transition could 
be continuous, meaning that the staging occurs progres- 
sively in the sample. It is indeed well known that GICs 
are extremely sensitive to these kind of transitions oc- 
curring isothermally under very modest pressures. For 
example KC24, which is a stage 2 GIC, starts a staging 
transition at 2.5 Kbar versus a stage 3 structure. The 
transition is continuous and completely achieved at ^ 6.5 
Kbar^^. Similar transitions have been reported in refi^ 
for MCg and MCi2n (n=2,3,4) with M=Rb, Cs. If such 
a transition occurs in CaCe then it would not necessarily 
show up in the Tc versus pressures curves, but it would 
explain the disagreement between experiment and the- 
ory. Of course other more complicate explanations are 
possible, including unconventional superconductivity. In 
all the case, the origin of superconductivity in interca- 
lated GICs is not yet completely understood. 
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